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Abstract: We present new quasi-exactly solvable models with inverse quartic, sextic, 
octic and decatic power potentials, respectively. We solve these models exactly by means 
of the functional Bethe ansatz method. For each case, we give closed-form solutions for the 
energies and the wave functions as well as analytical expressions for the allowed potential 
c~| \ parameters in terms of the roots of a set of algebraic equations. 
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1 Introduction 



> 

Since the early works on singular potentials (see e.g. pQ- [3]), an extensive literature has 
been developed on the subject. By singular potentials, we mean those potentials V(r) 
with property lim r _>o r 2 V(r) — > oo, although sometimes the inverse-square potential is 
| also regarded as singular. The investigation of singular potentials covers a wide range of 

physical and mathematical interest. In view of the availability of a comprehensive review 
article by Frank, Land and Spector [1] on singular potentials and applications, here we 
^ . will only refer to the main points of the topic. 

' One of the early works that generated much interest in the study of singular potentials 



was the one by Predazzi and Regge [5] , who argued that real world interactions were likely 
to be highly singular and thus the study of singular potentials rather than regular poten- 
tials might be more relevant physically. This was thereafter followed by the applications of 
the singular potentials r~ n (n > 2) in the study of the (p, p) and (p, ir) processes in high 
energy physics @ [7] . The interactions of nucleons with K-mesons and a-a scattering have 
been reproduced by repulsive singular potentials [5]. As examples from non-relativistic 
quantum mechanics, the problem of high-energy scattering by strongly singular potentials 
was investigated by many authors (see e.g. [H1QD, Hj)- In molecular physics, inter-atomic 
or intermolecular forces are mostly represented by singular potentials whose parameters 
are determined phenomenologically. Examples are the Mie-type potential [12t [13] and 
the Lennard- Jones potential r -12 p3] which describe the interaction between diatomic 
molecules and two non-polar molecules, respectively. 
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In field theory, the importance of singular potentials emerged from the efforts to find 
effective potentials of field-theoretic interactions in the Bethe-Salpeter equations. In [15] , 
Bastai, Bertocchi, Fubini, Furlan and Tonin discovered a remarkable relationship between 
the renormalizability of a field theory and the regularity of the effective potential. That 
is, the effective potentials for non-renormalizable field theories are singular, whereas su- 
perrenormalizable and renormalizable field theories give rive to regular and "transition" 
potentials, respectively. Thus, any new insight gained in the analysis of singular potentials 
could lead to a better understanding of quantum field theories which are not perturbatively 
renormalizable. 

In this paper, we consider a class of most frequently discussed singular potentials 
in nonrelativistic quantum mechanics, that is the spherically symmetric, inverse power 
potentials of the form [U [16j [T7] , 

k=0 

where r £ (0,oo), at and are positive real numbers. The corresponding Schrodinger 
equation for the radial wave function ^(r) is given by 



d 2 tU + 1) 22 . . 



#(r) = 2Etf(r), w > 0, 



where I = — 1,0,1, and E is the energy eigenvalue. Much of the investigation of 
singular potentials is concerned with the solutions of this equation. Several approximation 
techniques are available in the literature for calculating the eigenvalues and wave functions 
of some inverse-power potentials (see e.g. [HI [191 USB ED 122] [23]). 

We show that at least four classes of singular potentials in the above family are quasi- 
exactly solvable, i.e. have polynomial solutions |24 [ I25 [ [26]. Namely, we present new quasi- 
exactly solvable models with inverse quartic, sextic, octic and decatic power potentials, 
respectively. We solve these models exactly by using the functional Bethe ansatz method 
presented in [27] . For each model, we obtain closed-form solutions for the energies and 
wave functions as well as analytic expressions for the allowed potential parameters in 
terms of the roots of a set of algebraic (Bethe ansatz) equations. To our best knowledge, 
our results on the singular quartic and sextic power potentials are largely new, while our 
results on the inverse octic and decatic power potentials are completely new. 



2 Quasi-exactly solvable inverse quartic power potential 

The inverse quartic power potential 

h r d 

d>0 2.1 



rr , . a b c d 
V (r ) = - + - + - + -j, 



was investigated by Predazzi and Regge [5] to determine analytic properties of the scat- 
tering amplitude in the case of a singular potential. The corresponding radial Schrodinger 
equation is given by 



d 2 £U + 1) „ 9 2a 26 2c 2d 



*(r) = 2E^(r) 



(2.2) 



Phenomenologically, this type of singular potential is a very useful form of anharmonicity 
in physical applications [20J. 
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Here we are interested in finding exact solutions of the Schrodinger equation. To this 
end, we first extract the appropriate asymptotic behaviour of the wave function \P(r) by 
making a substitution. After a brief inspection of the differential equation, we arrive at 
the transformation, 



\fr(r) = r 7 exp 



9 D 

Ar 2 +Br-\ 

r 



x(r), 



(2.3) 



where 7, A, B and D are constant parameters to be determined. Substituting into the 
differential equation, we find that acceptable asymptotic behaviour of \I/(r) requires 



7=1 + 



> 0, A 



UJ 



D 



2d. 



'2d 2 

At this stage, the parameter B is free. Now the differential equation for x(r) reads 

* r + V2d 



(2.4) 



rV (r) + 2 



-wr + Br z + 1 + 



x'(r) 



+2 



B 1 + 



-!- { -2ojBr 6 + 
a — uj\[2d 



2E + B 2 - u 3 + 



2c 



r }> x(r) = (2b + £(£ + !)- AB\/2dj x(r) 



(2.5) 



In the following, we show that if B = or uj = 0, then Eq. (|2.5p is quasi-exactly solvable 
provided that the potential parameters satisfy certain constraints, and exact solutions are 
given by the degree n polynomials 



x(r) = Y\_( r ~ r *)' x ( r ) = 1 f° r n = 0- 



(2.6) 



8=1 



Case 1:5 = 

In this case, Eq. (|2.5|) becomes 



r 2 x"(r) + 2 



-wr J + 1 + 



2E - oj 3 + 



2c 



r 2 -2 



2d 

a + wV^rf 



r + v 7 ^ 



aj'(r)+ 



r}x(r) = [2b + £(£ + l)]x(r). (2.7) 



Substituting (j2.6|) into (j2.7|) and applying the functional Bethe ansatz method in [27] 
which is outlined in the Appendix, we obtain the closed-form expressions for the energies 
and the wave functions 

" -• 1 1 + \ + -k.. 



n« 



r — rs 



and the constraints for the potential parameters 

/ n 



U 2 

exp I --r 



2d 



(2.8) 



-uj 



(2.9) 



i=l 



3 



(2 \ n 
n+l + ^L) -2uVr? = 2ft + ^ + l), (2.10) 
V2d/ 



where the roots {r{\ satisfy the Bethe ansatz equations 



^ = * ^ , z = l,2,...,n. (2.11) 

The wave functions ^ n (r) are squarely integrable, i.e. J °° \^ n {r)\ 2 dr < oo. In fact the 
normalization constants of ^f n (r) can be calculated analytically by means of the standard 
integral 



2 M2\ , _ $ r,io (1 1 v v 



r v exp ( -mr - — )dr= ^_ G 02 



2 u + 1 ^¥jT[ uz \2' 2 2' 2 



1 2 



for Re(z/) > 0, Re(/ii) > and Re(^) > 0, where G™ n is the Meijer G-function. 

As examples of the above general expressions of the exact solutions, we study the 
ground and first excited states of the system in some detail. It is easy to see that x = 1 is 
a solution of (|2.7p for certain values of the potential parameters. This solution corresponds 
to the n = case in the general expressions above. Indeed, from Eqs. (|2,8p and (|2.9p . 
([240]) . we obtain 

f% c 
£ = u - + 



.2 V2d 

a = -ujVm and 2b = -£(£ + 1). (2.12) 
The corresponding wave function is 

* (r) = r 1+c /^exp ^-|r 2 - . (2.13) 

This wave function has no nodes and so the state described by it is the ground state of 
the system. 

The first excited state solution of the system corresponds to the n = 1 case of the 
general expressions above. Explicitly, the energy and wave function read, respectively, 

E l =uj[- + 



.2 y/2d, 

^(r) = r 1+c /^(r - r x ) exp I --r 2 - ^ | (2.14) 
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with the constraints for the potential parameters 

a = —uj ( 7*i + V2d) 



'2d. 

Here the root ri is determined by the Bethe ansatz equation 



wr? - ( 1 + -£= ) n - \/2d = 0, 
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giving rise to 



1 

2uj 



1 + 



1 + 



+ 4wV2d 



(2.16) 



Case 2: u = 

As A = in this case, in order for the wave function ^(r) to have acceptable asymptotic 
behaviour, we require the constant parameter B to be negative. 

The radial Schrodinger equation corresponding to this case was previously investigated 
in |28| (see also |29| for the ground state solution of the Schrodinger equation). Here we 
solve this system exactly by using the functional Bethe ansatz method of 

If we set 



2E + B 2 



0. 



(2.17) 



then Eq. (|2.5|) becomes 



+2 



r 2 x"(r) + 2 



B 1 + 



Br 2 + 1 + 



r + V2d 



x'{r) 



rx(r)= 2b + £(£ + 1) - 4BV2d x(r) 



(2.18) 



Substituting (|2.6p into (|2.18p and following the procedure in [27J (or applying the results 
of the Appendix), we obtain the relations, 



B 



n + l + 



2d 



2B 



2V2d + r ij + n(n + l + -jL^J =2b + 1{£ + 1) 



(2.19) 
(2.20) 



and the Bethe ansatz equations 



E 



+ 



Br 2 +(l + j=)r l + V2d 



,.n- rj 



1,2,..., n. 



(2.21) 



As mentioned above, we require B < 0. It follows from ()2.19p that the potential 
parameter a has to be negative. So we assume that a < if uj = 0. 

Thus we have the closed form expressions for the energies and the wave functions, 



E„ 



2 \ n + 1 + w d , 



.1=1 



cxp 



tf n (r) = r 1+c ^ 2d 
and the constraint for the potential parameters 

2a 



2d 



n + l + ir d 



^ + § ri ) +B ( B + 1 + |) 



n + l + 



(2.22) 
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n + 1 



[2b + £(£ + 1)] 



(2.23) 



with the roots {rj} determined by (|2.21|) . Again the wave functions ^ n (r) are squarely 
integrable, and the normalization constants can be computed analytically with the help 
of the standard integral [33] 



r exp y—^ir J dr = 2 I — 



BesselK [y + 1, 2y7*l7*2) 



for Re(z^) > 0, Re(/ii) > and Re(^ 2 ) > 0. 

Similar to the Case 1 above, the ground state energy and wave function are obtained 
from the general expressions by setting n = 0, 



1 



2 V 1 + -r=, 

\ V2d, 



En 



T i \ i+c/s/2d / ay 2d 
$ (r) = r i+c/v/d exp [ ^=r 



2d 



c + V2d 



with the constraint for 



8ad 



2d. 



(2.24) 



(2.25) 



2b + £{£ + !) 

The first excited state solution corresponds to the n = 1 case with the energy and the 
wave function given by 

/ \ 2 

1 / a \ 



2 \ 2 



2d . 



^i(r) = r 



i+c/v/a (r - n) exp I — r 



2d 



(2.26) 



x c + 2V2rf 

where the potential parameters satisfy the constraint 

2c 2 - cV2d [2b - 6 + £(£ + 1)] + Ad a (n + 2y/2d\ - 2(6 - 1) - l{£ + 1) =0 (2.27) 
and the root r\ is determined by 

9 

arf 



2 -i 



+ 1 + 



^ n + V2d = o 



c + 2-\/2d 



2aV2d 



1 + 



1 + 



8ad 



c + 2V2d 



(2.28) 



3 Quasi-exactly solvable inverse sextic power potential 

Pais and Wu |9] studied the problem of scattering by the singular potential d/r 2+2n + 
e/r 2+n (n > 1) in non-relativistic quantum mechanics. In this section we will consider the 
n = 2 case, i.e. the inverse sextic power potential 



(3.1) 
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This potential has been used in atomic, molecular and nuclear physics [16\ I17j . The 
corresponding radial Schrodinger equation is 



dr 2 



1(1+1) o 2 2e 2d 



tf(r) = 2E^(r) 



(3.2) 



The exact ground state solution of the Schrodinger equation was previously obtained in 
|31| . Here we review some of the results given in [30] and present the general exact solutions 
of the system. 

In order to find the exact solutions of the Schrodinger equation, similar to the inverse 
quartic power potential case we extract the appropriate asymptotic behaviour of the wave 
function Vl/(r) by making the substitution, 



^f r \ = r 3/2+e/V2d 



exp 



— r z 



2d 1 



v(r), 3/2 + e/V2d > 



(3.3) 



we obtain the differential for v(r), 



„, x 2 / 2 3 e V2d \ ., . 
v"(r) + - -ur 2 + - + -= + —5- v'{r) + 2 
r\ 2 ^2d r 2 ) w 



2ojV2A + (l + 



2rf 



J5-w 2 + 



+ 1 



u(r) 



u(r) 



(3.4) 



Then the change of variable t = r 2 transforms (|3.4p into the form, 



*V(i) + 



-ujt 2 + 2 + 



t + V2d 



At) + 



E-uj\2 + 



v(t) 



2uV2d + (l + 



2d 



+ 1 



v(t). 



This equation has the degree n polynomial solutions 

n 

v(t) = Y[(t - ti), v(t) = l for n = 



(3.5) 



(3.6) 



i=i 



with distinct roots t{ provided that the potential parameters satisfy certain constraints 
[30J. The closed- form expressions for the energies and the wave functions are 



E n = u>[2n + 2 + 



* n (r) = r 3 / 2 + e /^ M 



i> 2 - 



cxp 



UJ 2 

-— r 
2 



2d 1 



(3.7) 



and the constraint for the potential parameters reads 

2 



2u (^/2d + 2 J^ti j + ^ + 



4n ( n + 1 + ] + ( + 1 



2rf 



(3. 
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while the roots {ij} are determined by the Bethe ansatz equations, 



ujtj - 2 + 



2d 



i = 1, 2, • • • , n. 



(3.9) 



The wave functions ^ n [r) are squarely integrable, and the normalization constants can be 
evaluated by using the standard integral 



2 ^2 



r v exp ( — /Uir ^ ) dr = ( — 



/'2 



(f+l)/4 



BcssclK 



z/ + l 



2y^i7^2 



for Re(z/) > 0, Re(/ii) > and Re(^ 2 ) > 0. 

As examples of the above general expressions for the exact solutions, we study the 
ground and first excited states of the system. The n = case gives the ground state 
energy and wave function 

E = uj ( 2 + 



3/2+e/V2d 



exp 



UJ 



2d 1 
~2~^ 



with the potential parameters constrained by 



2uV2d + ft + 



2d 



+ 1 



(3.10) 



(3.11) 



The first excited state solution corresponds to the n = 1 case of the general expressions 
Eqs. (|3.3p and (|3,6p - (|3.9p . The energy and wave function are given respectively by 



E x = uj 4 + 



UJ 9 

— r z 
2 



y i{r)=r 3/2+e/V2d (r 2 -tx) exp 
where the root t± is determined by the Bethe ansatz equation, 

ujt\ - ( 2 + -J= )h- \f¥d = 

1 



2d 1 
2~^2 



(3.12) 



2d 



2 + 



+ 4o;V2(i 



(3.13) 



and the potential parameters obey the constraint, 



e / — 4e 
— - 2ujV2d + 5 + -= 
2d v 7 ^ 



1 

+ 2 



2 + 



+ 4wV2d 



(3.14) 



Here (|3.13[) has been used in deriving this equation from (|3.8p . 
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4 Quasi-exactly solvable inverse octic power potential 



In this section, we consider the inverse octic power potential 



abode f q h 
- + ^T + ^T + -r + ^ + 4r + 47 + ^r 



h>0. 



(4.1) 



This potential is an extension of the inverse quartic and sextic power potentials presented 
in the previous two sections. We will show that, in spite of the presence of additional 
coupling constants, this potential remains to be quasi-exactly solvable. 
The corresponding radial Schrodinger equation reads 



d 2 1(1+1) 9 o 2a 2b 2c 2d 2e 2f 2g 2h 



*(r) = 2EV(r). 

(4-2) 

To obtain exact solutions, we need to extract the asymptotic behaviour of the wave func- 
tion Vl/(r) and transform the differential equation into an appropriate form. Our experience 
with the the previous two simpler cases suggests the following substitution for ^(r): 



\P(r) = exp 



D F 



Ar + Br H h 



r r^ 



G 



y(r), 



(4.3) 



where f3, A, B, D, F and G are as yet unknown constant parameters. Substituting into 
the Schrodinger equation and after some algebras we find that the acceptable asymptotic 
behaviour of ^f(r) requires 



A 
F 



2 + ^ + 9 



g 



OJ 

2' 



y/2h 4 V h? \Ah 



/ >0, 



D 



9 



2V2h 



G 



2h 
1 



/ 



9_ 

Ah 



2h. 



(4.4) 



The parameter B is still free at this stage. Now the differential equation for y(r) is given 
by 



rY(r) + 2 



-ur 5 + Br 4 + /3r 3 + 



/ 



Ah 



+ 



gr 



+ <j -2uBr 5 + (2E + B 2 - u - 20u) r 4 - 2 

2B 



a- /3B + 



2h 

<jj 



+ y/2h 



y'(r) 



f 



Ah 



+ 



(p-l)(J3-£-l)-2b 
0-1 ' - 2 



2d 



2h 

9(2(3 - 3) 



f 9 2 \ 2guj 
2h V 4V J2h 



gB 



2h 



2BV2h 



1 

2h 



2h. 
f- 



r> y{r) 

o\ 2 

r_ 

Ah 



y(r). 



(4.5) 



If B = or oj = 0, then Eq. (|4.5|) is quasi-exactly solvable provided that the poten- 
tial parameters satisfy certain constraints, and exact solutions are given by the degree n 
polynomials 



yir) 



n< 

i=i 



r - n 



y(r) = 1 for n = 0. 



(4.6) 
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This is shown as follows. 



Case 1: B = 

In this case, equation (|4.5p becomes 



rV(r) + 2 



+ 



-cjr 5 + /3r 3 + 

+ <{ ( IE -u- 2(3uj) r 4 - 2 
2gw 



Co' 



2/i V 4/i 



(P-£)(P-£-l)-2b 



r 2 -2 



a + 



v 7 ^ V 4/i 



r y(r) 



g(2/3-3) 1 / g 
2/i V J 4/i 



2\ 2 



(4.7) 



Substituting (|4.6p into this equation and applying the results in the Appendix, we obtain 
the closed form expressions for the energies and the wave functions 



E n = LO 



5 e g 2 

n+ 2 + 7E + 4Vh?{4h- f 



Y[(r - n) 



i=l 



U 2 

exp I --T 



1 



g 2 \ 1 g 1 y/2h 1 



7 2/i V" 4/t/ r 2\/2hr 2 3 r 3 
with /3 given by (|4.4p and the constraints for the potential parameters, 

-,2 ' 



(4. 



a = —uj 



b=\[W + e)((3-e-l)+n(n + 2P-l)]-^L-u;^2rl 



c = -a,J>? + (n + /3 - 1) £r, + U1V2H+ (f - 

8=1 j=l 

i=l j=l i<j i=l 



(£_ 

4/i 



q(2n + 2/3-3) 1 / g 2 \ 2 
2^ + 4/il J 4/i J ' 



(4.J 



where the roots {r{\ satisfy the Bethe ansatz equations 



E 



+ 



0, i = 1,2, - •• ,n. (4.10) 



It can be checked numerically using Maple that the wave functions 9 n (r) are squarely 
integrable: / °° |\l/ n (r)| 2 dr < oo. 
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As examples of the general expressions for the exact solutions, we consider the ground 
and first excited states of the system in some detail. It is easy to see that y = 1 is a 
solution of (|4.7p provided the potential parameters satisfy some constraints. This solution 
corresponds to the n = case in the general expressions above. Indeed, from Eqs. (|4.3p . 
(03D, (USD, and we obtain 



-En 



2/i 4V/i 3 V 4/l 



t / \ fl 1^2 1 

Wo r = exp r . , . , 



(4.11) 



with the constraints 



b = ±[(13 + £)(?- 1- 1)]--jL, 



V2h V 4/1 



4h\ J 4/J + 2^ ■ 



(4.12) 



The wave function (|4.1ip has no nodes and so the state described by it is the ground state 
of the system. 

The first excited state solution of the system corresponds to the n = 1 case. From the 
general expressions above, we obtain the first excited state energy and wave function, 



2 V2h 4\h 3 \Ah 



*$>±(r) = r l:l (r — r\) exp I — — r 2 



1 



g 1 V2h 1 



>2h V" 4hJ r 2y/2h~r 2 3 r 3 
where the root r\ is determined by the Bethe ansatz equation, 

,2 ' 



(4.13) 



-url + Prf + ^=(f- fr-^l r\ + -%= ri + V2h = (4.14) 
\/2h V 4 V v 7 ^ 



and the potential parameters are constrained by 



-co 



£)((3 + £ + !)] 



ujg 



ujr 



2h 



l) 



P -(f- 9 -) 



-urf + j3r\ - wv2/i 



4/; 



+ 



5(2/3 - 1) 



2V2h 



(4.15) 
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Case 2: u = 



In this case, A = —uj/2 = 0. Thus in order for the wave function \P(r) to have 
acceptable asymptotic behaviour at r — > oo, we require the constant parameter B to be 
negative. Furthermore, we set 



2E + B Z = 0. 



(4.16) 



Then equation (|4.5p becomes 



r 4 y"(r) + 2 



5r 4 + /3r 3 + 



4V 



y'(r) 



+ ^ 2 - a] r 3 + 



(/3-^)(/3-£-l)-26 + 



2B 



J Ah 



2h 



Ah 



gB 



2d-&ZV-2BV2h 



2h 



2h. 
1 

2h 



r } y(r) 

2\ 2' 



/ 



EL 

Ah 



yir) 



(4.17) 



Substituting (|4.6|) into this equation and applying the results in the results of the Ap- 
pendix, we obtain constraints for the potential parameters, 



B 



n + /3' 



b= - [(P + £)(P-£-l) + n(n + 2(3- 1)] + B 



c = B 



+E ; 

i=l 



+ (n + /3-l) 



n 



Ah) 



i=l 

Ah V 4/i / 2^ 



n n n i/2\ n 

d = 5^rf + (n + /3-l)^r J 2 + ^r J r, + - 7 =( / -|-)^ 

i=l i = l i<j V \ / 



(4.18) 



where /? is given by Eq. (|4,4p and the roots {r^} satisfy the Bethe ansatz equations 



£ 



i , Br t + ^ + 7k{f-h) + ^ h + V2h 



+ 



0, n= 1,2,- •• ,n. (4.19) 



It follows that the energies and the wave functions of the system are 

2 



* n (r) = r p 



Y[(r - n) 



i=l 



exp 



1 



1 



2 \n + /3J ' 



5 1 v 7 ^ 1 



n + /3 y 7 ^ V 4/i/ r 2V2h~r 2 3 r 3 



(4.20) 
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As mentioned above, B is require to be negative. To satisfy this requirement we see 
from (|4.18p that the potential parameter a has to be negative. So we assume that a < 
if uj = 0. Then it can be checked numerically that the wave functions ^ n (r) are squarely 
integrable. 

Similar to the Case 1 above, the ground state solution is obtained from the general 
expressions by letting n = 0, yielding 



En 



a 



with the constraints 



r exp — r 
\ P 



i- 3 -Y- 

J Ah r 



9 



1 



2V2h~r 2 



2h 1 

3 r 3 



(4.21) 



h/3(p-i)-£(e + i)} + —?= 



J Ah 



ag 



03-1) 



d 



/V2/i + y/2h 

1 / 2 \ 2 

4/i V Ah 



f 



Ah 



aV2h g{2j3 - 3) 



P 2^2h 

The first excited state solution corresponds to the n = 1 case. We thus have 

2 



Eh 



l + P 



(4.22) 



^i(r) = r^(r — ri) exp 



subject to the constraints 



1 + /3 



J 4/i/ r 



5 



2y/2h~r 2 



2h 1 

3 r 3 



(4.23) 



El 

Ah 



d 



1 + 



i 



P4 + 



l + P 



rj + 0n 



ag 



(l + (3)V2h 



+ 



P 



Ah 



9 2 \ 



1 



/-7T +ir /-TT 



4/i7 4ft 



4ft 



aV2ft 5 (2/3 - 1) 



/9 + 1 



Here the root r\ obeys 
a 



l + P 



r\ + Pr{ + 



El 

Ah 



r\ + -4=n + V2h 



2h 



2y/2h 



0. 



(4.24) 



(4.25) 



Let us remark that the results obtained in this section are reducible to those of the 
inverse quartic potential when g, f, e, h — > so that (3 — > 1 + c/y2d and -j= If — |r 



/2/i v 4 V 

v2d, and reproduce those of the inverse sextic potential if a, b, c,e,g,h — > such that 
P -»• 1 + d/V27 and g/v 7 ^ -> v 7 ^/- 
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5 Quasi-exactly solvable inverse decatic power potential 

We consider the inverse decatic power potential 



_ a b c d 

V\ r ) — ^4 + ^6 + ^8 + ^10 ' 
where d > 0. The corresponding radial Schrodinger equation is given by 



d 2 £(£ + 1) on 2a 2b 2c 2d 



(5.1) 



(5.2) 



Similar to the cases studied in the previous sections, finding exact solutions to this equation 
requires first extracting appropriate asymptotic behaviour of the wave function \I>(r). This 
is achieved by means of suitable transformations of ^(r). After some algebras we arrive 
at the transformation 



\]/(r) = r v exp 



UJ 



c 1 V2d 1 



2V2d~r 2 4 r 4 



u(r), 



5 b c 2 12 



+ 



+ 



' 2 ^2d 16 V d 3 
Then the differential equation for u(r) is 



u"(r) + 



> 0. 



2d c ,i ,, , 

H — - — == + f] — ojr J u \t) 



(5.3) 



r 2 V2d 



+ 



77(77 - 1) - ^ + 1) - 2wc/V2d 3c + 2<4j - 2cr/ + 2aV2d 



r z V2(ir 4 
+2£u(r) = 0. 

A further variable change z = r 2 transforms the above equation into the form 



u(r) 



(5.4) 



z 6 u"{z) + 



-ujz i + (77 + 1/2) z 2 + -£=z + v 7 ^ 



u'(z) + |i[S- w (i / + l/2)]z 2 



77(77 - 1) 1) 



2wc 



1i(z) = - 

v ; 2 



a + u)V2d + -^=(3/2 - 77) 
v2d 



(5.5) 



This equation is quasi-exactly solvable and has the degree n polynomial solutions 

n 

s) = JJ(z - Zi), u(z) = 1 for n = (5.6) 



1=1 



with distinct roots Zj provided that the potential parameters satisfy certain constraints. 
Indeed, substituting (|5.6p into (|5.5p and applying the functional Bethe ansatz method in 
|27] outlined in the Appendix, we get the closed form expressions for the energies and the 
wave functions 

E n = .(2n + ri+ l)=uL + S + ^ d + f 6] [I), 
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»„(r) = r" 



n< 

.i=i 



r 2 - Z,; 



Co' 



6XP "2 r "^5r» 



1 \/2d 1 



4 r 4 



(5.7) 



and the constraints for the potential parameters 

\ 2 



2co\ -4= + 2j>) = (^-^l + 2] _^ + IV+4n(n + 2+ C ' 



2d 



16 V d 3 



a = -2a;^z l 2 + ^4n + 2 + ^Aj^ Zi + _^ (2n + i) + 
where the roots {z{} obey the Bethe ansatz equations 



wV2d, (5.8) 



j -wv^d ^ 3 + (3>/2d + c 2 /8d) zf + czi + 
+ f=-^ 



2d 



2d zf 



= 0, i = 1,2,- •• ,n. (5.9) 



The wave functions ^ n (r) are squarely integrable, i.e. J °° (^(r^pdr < 00. This can be 
checked numerically using Maple. 

Clearly, the ground state solution can be obtained from the above general expressions 
by letting n = 0, yielding 

fn b C 

Eq = u\3 + + 



^o(r) = r 



5 , 6 



2" r V2d + 16 V 



v 7 ^ 16 V d 3 ^ ' 
/ u 2 c 1 \/2d 1 ~ 



(5.10) 



with the constraints 



2c 



2d 



16Vd3 +2 



4 



c c 

+ 



6jV2d 



The first excited state solution corresponds to the n = 1 case. We thus have the energy 
and the wave function, 

b 



E 1 =u\5 + 



c 2 12 



5 1 6 c 



vp^r) = r 2 ls V?(r 2 — 



V2R ie\ d 3 )' 

c 1 v 7 ^ 1 



N ^ 2 

zi) exp — — r — 



, 2' 2V2dr 2 4 r 4 ^ 
with the potential parameters subject to the constraints 



uj (c + 2v / 2dzi) 



16 V d 3 +2 



'4 



a = -2ujzj + 6 + 



c 2 /2 



3c c 3 

«1 + -F= + 



8 V d?) 1 ^2d 16d 2 
Here the root z\ is determined by the Bethe ansatz equation, 

-uV2d zf + (3\/2d + c 2 /8d) z 2 + cz\ + 2d = 



+ 6\/2d + ^, 
4d 



- w\/2d. 



(5.12) 



(5.13) 



(5.14) 
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6 Concluding remarks 



We have presented four classes of quasi-exactly solvable anharmonic singular potentials 
and solved the corresponding radial Schrodinger equations exactly by means of the Bethe 
ansatz method. Closed-form expressions for the energies and the wave functions of the 
singular models have been obtained for certain values of the potential parameters. Analytic 
expressions for the allowed potential parameters are given explicitly for each of the four 
cases in terms of the roots of a set of algebraic Bethe ansatz equations. 

It is remarkable that the four classes of anharmonic singular potentials are quasi- 
exactly solvable. We have thus added four new members to the club of quasi-exactly 
solvable quantum mechanical systems. It is hoped that our findings in the present paper 
will lead to new applications for these singular potentials. 
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A Appendix: the functional Bethe ansatz method 

As seen from (|4.7p and (|4.17p . the inverse octic power potential leads to the 2nd order 
differential equations of the form, 

S(t) = 0, (A.l) 
where P(t),Q(t) and W(t) are polynomials of at most degree 4, 5 and 4 respectively, 

4 5 4 

P(t) = J2Pkt k , Q{t) = Y.1^^ W(t) = J2^kt k , (A.2) 

fc=0 k=0 k=0 

Pk,Qk an d Wk are constants. This 2nd order differential equation is a generalization of 
the one considered in [27] and does not seem to have ever appeared previously in the 
literature. A similar differential equation also appears in quantum mechanical models 
with (non-singular) octic potentials [32] . 

In this appendix, we show that the differential equation (|A.1|) is quasi-exactly solvable 
for certain values of its parameters and exact solutions are given by polynomials of degree 
n in t with n being non-negative integers. We apply the functional Bethe ansatz method 
in |27j . We seek polynomial solutions of the form 

n 

S(t) = Y[(t-ti), S(t) = l for n = (A.3) 

i=l 
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with tx,t 2 , . . . , t n distinct roots to be determined, then Eq. (|A.lj) becomes 

n 1 n 2 
-w = [pit 4 +p 3 t 3 +p 2 t 2 + Pi* +Po) Et~tE 



t=l 

n 



+ {Q5t 5 + q±t A + q 3 t 3 + q 2 t 2 + qi t + q ) V 7^- 

i=l f Ji 

+w 4 i 4 + W3^ 3 + + Wit (A.4) 



The left hand side of this equation is a constant, while the right hand side is a meromorphic 
function with simple poles t = i, and singularity at t = 00. For this equation to be valid, 
the right hand side must also be a constant. We thus demand that the coefficients of 
the powers of t as well as the residues at the simple poles of the right hand side be zero. 
Following Liouville's theorem, this is the necessary and sufficient condition for the right 
hand side of (|A.4j) to be a constant. It follows that 

(A.5) 
(A.6) 



(A.7) 



W4 = 


-nQ5, 

n 






U> 3 = 


~95 E " 

i=l 






U>2 = 


n 

"95 J>i 

i=l 


n 

- 94 E ** - 

i=l 


n(n - l)p 4 - n<? 3 , 


W\ = 


n 

"95 

i=l 


n 

-94E*i " 

i=l 


n 

[2(n - l)p 4 + q 3 ] y~]U- n{n - l)p 3 
i=i 


W = 


n 

"95 E^ 

8=1 


n 

-94E*i " 
i=l 


+ 2(n - l)p 4 ] E *i " 2 ^ 4 E 

8=1 i<J 



- [2(n - ljpa + 92] E^ " n ( n " 1 )^ 2 " n 9i' ( A - 9 ) 
i=i 

where the roots ti, t 2 , ■ ■ ■ , t n satisfy the Bethe ansatz equations 

V 2 + 95tf + g4 + g3*| + gag + qiU + 9o = Q 2 r f A _ 

^.U-tj Pitf + P3 t 3 + p 2 t 2 + Pl ti + po ' 

Equations ()A.5p - (jA.9p give all polynomials W(t) such that (jA.ip has degree n polynomial 
solutions (|A.3p , The above method can be easily generalized to differential equations with 
higher degree polynomials P(t), Q(t) and W(t), provided that degW(t) < deg Q(t). 
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